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Abstract 

We considered the thermodynamics in spaces with deformed com- 
mutation relations leading to the existence of minimal length. We 
developed a classical method of the partition function evaluation. We 
calculated the partition function and heat capacity for ideal gas and 
harmonic oscillators using this method. The obtained results are in 
good agreement with the exact quantum ones. We also showed that 
the minimal length introduction reduces degrees of freedom of an ar- 
bitrary system in the high temperature limit significantly. 
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1 Introduction 



Now it is widely accepted that quantum gravity and string theory lead 
to minimal length, i.e. there must exist a lower bound to the possible 
resolution of the distance [1-8]. The minimal length appears due to a 
modification of the standard Heisenberg uncertainty relation: 

MAP > ^(1 + (3AP 2 + ...), (1) 
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where (3 is a positive parameter. This modified relation implies that 
uncertainty of coordinate AX is always larger than AX m i n = h\/~/3. It 
was shown that such a modified uncertainty relation can be obtained 
as a natural consequence of a specific modification of the usual com- 
mutation relation [9-13]. We would like to note that modification of 
the commutation relation is only one of the possible ways of introduc- 
ing minimal length. For a review of different approaches to theories 
with the minimal length see [14] . 

The deformed commutation relation according to Kempf et al. [11- 
13] may read in one-dimensional space 

[X,P]=zh(l + (3P 2 ) (2) 

or in D-dimensional space (i, j £ 1 . . . D) 

[X t , Pj] = ih{(l + pP 2 ^ + tfPiPj), [Pi, Pj] = 0, (3) 

l X i ,X l] ^-^l 2 ^ P \ p i X j - P} X i) . (4) 

It is believed widely that the minimal length can be helpful in 
the context of divergencies regularization in different areas of theoret- 
ical physics [5]. Contrary to the previous theories, where the minimal 
length appears phenomenologically, in our case the minimal length ap- 
pears as a natural consequence of modified commutation relations. An 
explicit representation of position and momentum operators is known 
and the existence of possible regularization can be verified from the 
first principles. It was shown that deformed commutation relations 
really might lead to regularization in quantum field theory [15]. Regu- 
larization with the help of deformed commutation relations for specific 
eigenvalue problems was studied in [16, 17]. 

Other implications of non-zero minimal length were considered 
in the context of the following problems: harmonic oscillator [12, 
13, 18-21], hydrogen atom [22-25], gravitational quantum well [26], 
the Casimir effect [27, 28], particles scattering [29, 30]. In several 
of the mentioned papers attempts to estimate the value of minimal 
length upper bound were made by comparing theoretical predictions 
and experimental data. In [22-25] considering the Lamb shift the 
authors estimated Al m j n < 10~ 16 . . . 10~ 17 m, analysis of electron 
motion in a Penning trap also gives AJ min < l(T 16 m [18], consid- 
eration of neutron motion in the gravitational field [26] gives large 
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AX m i n < 2.4 • 10 _9 m due to significant experimental errors. In our 
analysis we will use AX m i n = 10 -16 m. 

Consideration of some theoretical physics topics requires the use of 
statistical methods. The purpose of this paper is to develop statisti- 
cal physics with some general form of deformed commutation relation. 
There are several papers on the subject. In [31] the statistical physics 
is constructed by modification of elementary cell of space volume ac- 
cording to modification of the commutation relation. In that paper a 
very strong restriction [JQ,-X,-] = was assumed. A similar analysis 
was performed in [32,33]. Density of states and black-body radiation 
were considered in [34] with the help of Liouville theorem analog. The 
microcanonical approach to ideal gas was considered in [35]. 

In this paper we develop a method for the consideration of ther- 
modynamical properties of the system with arbitrary commutation 
relation between momenta and positions operators. We approximate 
the value of the partition function with the help of a semiclassical ap- 
proach. We may expect that such an approximation is applicable since 
the WKB approximation gives correct predictions for one-dimensional 
eigenvalues problems with deformed commutation relations [36,37]. 
We apply this method to the particular case, namely to the three- 
dimensional Kempf's commutation relation ([3H1]) and consider impli- 
cations of the minimal length on thermodynamics. 



We consider a system consisting of N identical non-interacting par- 
ticles in the external field U. The behavior of each particle may be 
described by Schrodinger Hamiltonian 



We use this Hamiltonian to analyze examples, but the developed 
method can be applied to a system described by Hamiltonian of any 
form. 

In this paper we assume that the partition function of a quantum 
system can be evaluated as 



where E n are eigenvalues of the Hamiltonian, n denotes different 
states. This formula has the same form either in the deformed or 
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(5) 




(6) 



n 
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in the non-deformed cases. In the non-deformed case one can use a 
semiclassical approximation for the partition function 

Z = J e- H ' T {dx)(dp), (7) 

where (dx) means dx\dx2 ■ ■ ■ dxp- Since it is the non-deformed case, 
the coordinate and the momentum variables Xi and pj are canonically 
conjugated {xi,pj} = Sij, {xi,Xj} = {pi,Pj} = 0. Here and below 
small variables x, p denote canonically conjugated variables. 

We show below that the existence of deformation becomes per- 
ceptible only for very high temperature. For high temperature the 
difference between Bose-Einstein, Fermi-Dirac or Boltzmann statis- 
tics is irrelevant, therefore we use Boltzmann one as the simplest. 
One of the aims of this paper is to generalize formula (|7|), which gives 
correct predictions for high temperature, to the case of the deformed 
commutation relation. Peculiarities of Bose-Einstein and Fermi-Dirac 
statistics for the particular form of deformed commutation relations 
were considered in [32,33]. 

Let us consider the general case of deformed commutation relation 

[X i ,P j ]=ihf ij (X,P), (8) 

[P i ,P j ] = ihh ij (X,P), (9) 

[X i ,X j ]=ihg ij (X,P). (10) 

In the above expressions operators Xi, Pj correspond to the same par- 
ticle. Operators describing different particles commute. The deforma- 
tion functions fij, gij, h, L j are restricted according to the properties 
of commutators: bilinearity, the Leibniz rules and the Jacobi iden- 
tity. The investigation on this subject for some special cases of the 
deformation function can be found in [13, 15]. 

In the classical limit fi ^ the deformed commutation relations 
flglfTOD lead to the deformed Poisson bracket [34, 38, 39] 

{X u Pj} = fij(X,P), {Pi,P 3 } = hij(X,P), {Xi,Xj} = 9ij (X,P).(ll) 

These Poisson brackets possess the same properties as the quantum 
mechanical commutators (IBT llOp . namely, they are anti-symmetric, bi- 
linear, and satisfy the Leibniz rules and the Jacobi identity [38,39]. 
According to Darboux theorem [40] it is always possible to choose 
such auxiliary canonically conjugated variables Xi and pi that Xj and 
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Pi as functions of x\ and pi satisfy equations (|TT|) . Where, just as an 
example, 

ftf;1 = ^5^« = fc( , P) . (12) 

^ dx k dp k dx k dp k 

Then, the initial Hamiltonian © can be considered as the func- 
tion of (x,p) and the one particle partition function can be evaluated 
according to formula ([7]). In appendix [A] we show that the Jaco- 
bian J = can always be expressed as a combination of Poisson 

brackets (fTTT) . It is an important result since it gives a possibility to 
calculate the partition function without introducing canonically con- 
jugated auxiliary variables. Namely, 

c -h(x,p)/t jdX){dP) 

This equation is an analog of formula ([7]) for the deformed case. It is 
an essential and expected result that the partition function does not 
depend on auxiliary variables since J depends only on the structure 
of the Poisson brackets (jlip which, together with the Hamiltonian, 
define all the properties of the system. 

Expression (|13p can be considered as a semiclassical approxima- 
tion of partition function ©. It is well known that the semiclassical 
approximation is good for large quantum numbers (it corresponds to 
high temperature regime). We will observe this situation with the 
help of the harmonic oscillators ensemble example: heat capacities 
calculated according to formulae © and (fT3|) coincide for high T. 



3 Examples 



In this section we analyze the effect of deformation ([SHU) on the ther- 
modynamical quantities. Since Jacobian (|38|) does not depend on 
coordinates and Hamiltonian has the form ([5]) evaluation of the one- 
particle partition function f|13f) can be separated in two parts: inte- 
grating over (X) and over (P). Namely, 



Z= /(dX)exp 



U(X) 



(dp) 



cxp 



P 2 



2mT 



(l + (3P 2 ) 2 (l + ((3 + (3')P 2 )' 



(14) 



Note, that during the transition from quantum expression for the par- 
tition function ([6|) to the classical one (I14p . we implicitly assume that 
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/3 and f3' are independent of h, i.e. we keep the deformation param- 
eters fixed as h — > 0. Some comments on this topic can be found 
in [38]. 

It is easy to see that deformation changes the partition function 
perceptibly if [5mT ~ 1, P'mT ~ 1. For electron gas with the minimal 
length 10 -16 m (3mT ~ 1 if T ~ 10 18 K. The very high value of this 
temperature justifies our choice of Boltzmann statistics. 

For low temperatures (5mT <C 1, (3'mT <C 1 the partition function 
of one particle (|14p simplifies to 

Z = Z (l- 3(3/3 + P')mT + o(T)), (15) 

where Z = J (dX) exp [-U/T] J (dP) exp -P 2 /2mT is the parti- 
tion function of one particle described by the same Hamiltonian for 
the non-deformed case. The existence of deformation leads to a change 
of internal energy and heat capacity in such a way 

E = E - 3(3/3 + (3')NmT 2 + o(T 2 ), (16) 
C = C -6(3/3 + f3')NmT + o{T). (17) 

For large T: (3mT S> 1 the partition function of one particle, 
internal energy and heat capacity are 

E = E -^NT + O(l), C = C -^N + o(~\. (19) 

This result will be verified on the quantum level for two examples 
(ideal gas and ensemble of 3D harmonic oscillators). We see that 
the presence of deformation freezes three degrees of freedom at high 
temperature. In section 2] we show that such a freezing is a natural 
consequence of the minimal length introduction. 

3.1 Ideal gas 

Hamiltonian of one particle is H = P 2 . We consider that the entire 
system is confined in volume V, i.e. for each particle / (dX) = V. 
For low temperatures, the one-particle partition function reads 

Z = V{2nmT) 3/2 (l - 3(3/3 + f3')mT + o(T)). (20) 



6 




It leads to the following expressions for internal energy and heat ca- 
pacity 

E = ^NT - 3(3/3 + (i')NmT 2 + o(T 2 ), C = -N - 6(3/3 + (3')NmT + o(T). (21) 
In limit T — > oo the one-particle partition function is 

z _ V^ 2 A 1 ,_ 

^(^+VPTW) 2 \ 2rnT[3 ^]3TW " 

(22) 

Here we take into account the second term of the series expansion to 
obtain the second correction to the internal energy expression (the 
leading term, which must be proportional to T, is zero): 

E =2^ VF+W +o(1) - (23) 

The last equation can be explained with the help of the quantum anal- 
ysis. Although the problem of the particle in the box in the deformed 
case has not been solved exactly yet, the preliminary analysis shows 
that there exists only a finite amount of bound states [41,42]. For high 
temperature it must lead to the situation that all particles occupy the 
highest energy level and thus, heat capacity must tend to as we see 
from (1231) . 

It is interesting to note that equation of state is the same as in the 
non-deformed case and its form does not depend on the temperature 
value. Namely, it reads 

pV = NT. (24) 

3.2 3D harmonic oscillators 

p2 

One-particle harmonic oscillator Hamiltonian is, as usual, H = h 

2m 
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? X . For low temperature expression (|14p gives the following for- 
mula for the partition function 

Z= (^f) 3 (1-3(3/9 + p')mT + o{T)) (25) 

The eigenvalue problem of the harmonic oscillator has been solved 
exactly in [18] for the first time and the spectrum reads 



E nl = huj{ [n+-)\jl + m 2 uj 2 h 2 



m + 1) + m+£T 
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mu>h 



• (26) 



We calculate the one-particle partition function using formula ([6]) in 
linear approximation over /3, (3' exactly. This expression is rather 
cumbersome, but in limit h — > we obtain simple expression for it 

Z = ^-(l-3{3f3 + f3')mT) + 0(l/fi 2 ). (27) 

The leading term of the last expression differs from (|25p only by con- 
stant factor (2-Kfhf which does not have any effect on thermodynamics 
and which we omit in the initial definitions of approximation (J7J [T3j) . 
For large T expression (|14p gives 



z =(^y 2 7mrvFTw ( 1+0 (?))- (28) 

E= -NT + 0(1). (29) 

This result can be easily explained and verified as follows. The energies 
with large n give the main contribution to the partition function at 
high temperatures. In the deformed case for large n according to 

expression flUD E nl (f3 + f3')n 2 . Such an energy dependence 

on the main quantum number is the same as for a particle in the box 
for the non-deformed case. Thus, it is an expected result to reproduce 
particles-in-the-box internal energy dependence on the temperature 

Heat capacity dependence on temperature is plotted on Fig. [TJ For 
high T heat capacity tends to 3/2 what follows from (I29j) . One can 
see that approximation (fT3|) gives predictions very closely to the exact 
result for high T. For selected values of parameters (3 = (3' = 0.01, 
fi = u = 2m = 1 it was expected that effects implied by the deformed 

commutation relation become significant for T = — — ~ 100. As one 

pm 

can see this implication appears to be important for a much lower 
temperature, namely for T = 1 -j- 2. 
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4 Heat capacity at high temperature 
and the minimal length 

For a wide class of physically important systems, Hamiltonian can be 
expressed as a sum of kinetic and potential energy. In this section 
we estimate the contribution of kinetic energy to heat capacity in the 
limit T — > oo. 

Let us consider the non-deformed case. In the case of high temper- 
ature only large values of momenta contribute to the partition func- 
tion. For large p the Hamiltonian can be approximated as H = ap n 
(for Schrodinger Hamiltonian a = n = 2; for Dirac Hamiltonian 
a = c, n = 1). So, 

f°° <*p n n , D D f°° „ D-l , , 

Z(T^oo)~/ e~—p D - 1 dp = T- a -- e~ x X — dx (30) 

Such a dependence on T means that C(T — ► oo) = — , which repro- 

n 

duces several well-known results. 

Let us consider deformed commutation relations leading to the 
minimal length. In a one-dimensional case 

[x,P]=ih(i + f(x,P)) 

and Schrodinger uncertainty relation reads 

(/(X,P)) \ ft/1 (f(X,P)) \ 
AX -2{AP + AP ) -2{P + P ) ' 

where P = \J (P 2 ) > AP. Therefore, a nonzero minimal length exists 
if lim f(X, P) grows as P or faster for large momentum values. Then 

the Jacobian J = 1 + / for large P must grow as P or faster. In the 
D-dimensional case it must grow as P D or faster. 

Let us first consider the case J ~ P D and denote by Po such 
an intermediate value of momentum that H ~ aP n , J ~ ~/P D for 
P > Po, but H <t: T for P < Po. Then the partition function reads 

r, 1 H/T PD ~ ldP J P D ~ l dP 1 7 aP"/T dP 1 , T 

Z = / e~ H/T ; w / ; + - / e~ aP /T — = const+-ln — 

J J J J 7 J P 7 «P ( 

P 

(31) 

Heat capacity can be easily calculated and C(T — > cxd) = 0. 
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A similar procedure can be performed when the Jacobian grows 
faster than P D . It gives Z(T — > oo) = const and consequently C(T — > 
oo) = 0. So, in any case leading to the minimal length C(T — > oo) = 0. 

We may conclude that for a wide class of Hamiltonians being pro- 
portional to a power of momentum for large momentum values and 
deformed commutation relations leading to the nonzero position un- 
certainty in the limit T — > oo the partition function of the system is as 
InT (the Jacobian grows as P D ) or does not depend on temperature 
(the Jacobian grows faster). Case of the constant partition function 
was discussed in the context of the string theory from the duality 
point of view [43]. A similar effect of degrees of freedom reduction 
was observed for two model Hamiltonians with particular deformed 
commutation relations in the context of grand canonical approach [31]. 

5 Conclusions 

In this paper we developed a general method for the consideration of 
thermodynamical properties of an arbitrary system of non-interacting 
particles with deformed commutation relations. To calculate a clas- 
sical expression for the partition function we express the Jacobian of 
making a change from arbitrary canonically conjugated variables to 
the initial variables as the combination of deformed Poisson brackets. 

We calculate the partition function for two model systems, namely 
ideal gas and ensemble of harmonic oscillators and considered corre- 
sponding thermodynamical properties. We showed that such a derived 
partition function is in qualitative agreement with the exact quantum 
expression in the limit of high temperature. Quite interesting conse- 
quence of the minimal length introduction is decreasing of the heat 
capacities of these model systems to and 3/2 for high temperatures. 
We managed to show that such a reduction of the heat capacity oc- 
curs (when compared with the corresponding non-deformed system) 
for each system with arbitrary deformed commutation relations lead- 
ing to the minimal length. Namely, the minimal length introduction 
completely removes translation degrees of freedom for very high tem- 
perature. 
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A Expressing Jacobian as a combina- 
tion of the Poisson brackets 

Let us denote Xj = A21-1, Pi = A21, Aj derivative with respect to x; L 
we denote Aj^i-i, with respect to pi as Aj^i- Then 

D 

{Ai,Aj} = y](Aj,2k-iAj,2k — A i: 2kAj t 2k-i)- 
fc=i 

Let us prove the following identity 

d(X 1 ,P 1 ,...,X D ,P D ) _ 1 ^ 

d(x 1 ,p 1 ,...,x D ,p D ) ~2°D\. ^ £ n... l2D iA 1 ,A 2 |...iA 2D _ 1 ,A 20 |, 

(32) 

where e^... j 2D is the Levi-Civita symbol. The right-hand of this iden- 
tity is equal 

2D ^ 2D 



■-•Dpi £ u— « 2 £>{^-u' ^4i 2 } " • • • " {Ai 2D - 1 i A{ 2D } — 2Djj\ e n».*2D 

' ii,...i 2 D=l ' ii,...i 2 zj=l 

D D 
^ (Al,2ji-lA 2 ,2ji — Ai 1> 2j 1 Ai 2) 2j 1 -l) ■■■ ^ (^i 2 D-i,2j' D -lA 2D ,2j D ~ ^i 2D _i,2j D A 2D ,2j D 
Jl=l JD=1 

= £)T 51 £ il---i2D J ^h,2ji-lAi 2t 2j 1 ■ ■ ■ Ai 2D _ lt 2j D -lAi 2Dj 2j D - 

In the latter equality we take into account that the Levi-Civita symbol 
is antisymmetric with respect to any indexes permutation, thus the fol- 
lowing property holds ^ £ii...i 2 D^il,2ji-l>W.h = - 51 e U...i2D^U,2il^i2,2il-l- 

From the fact that 



£ h...i 2D Ai 1: 2j 1 -lAi 2i 2j 1 • • • ^i 2D _i,2.7' D -l^.i2D,2.7£> — 



ll,...l 2 D 
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det 



2J1-1 



^2,2ji-l -^-2,2ji 



\ ^2D,2ji-l ^-2D,2ji 



^jjj-l 



A, 



,2jo 



A2D,2j D -l A2D,2j D J 



and using the determinant properties it becomes obvious that in the 
last formula of equation (|33[) only terms with different j give the con- 
tribution to the final result. All the terms with ji ^ j 2 ^ ■ ■ ■ 32D 
are equal. The total number of these terms equals D\. Thus, 



D\ . 5 . ? 

Ji,— 3d »i,—*2r> 



E i 1 ...i 2D Ai 1 ,2j 1 -lAi 2t 2j 1 ■ ■ ■ A-i 2D _ lj 2j D -lAi 2Dt 2j D 



12D 



.i,2D-iA aD ,2D = det (Ay), (34) 



ll,...l2D 



which proves identity (|32p . 

The right-hand side of expression (|32p contains (2D)! terms, each 
of them is a product of D Poisson brackets. Due to the skew-symmetry 
of Poisson bracket some of the terms are equal and the total amount 
of terms can be reduced to (2D — 1)!! terms. Below we enlist Jacobians 
for one-, two- and three-dimensional cases. 

In a one-dimensional case the following expression 



d(X,P) 
d(x,p) 



{*,P} 



(35) 



is obvious. 

In a two-dimensional case 



d(X 1 ,P 1 ,X 2 ,P 2 ) 
d(x 1 ,px,x 2 ,p 2 ) 



{ A X , P X } { A 2 , P 2 } - { A X , A 2 } {P X , P 2 } -{ A X , P 2 } { A 2 , P X } 

(36) 



can be checked by hand. 
A three-dimensional case: 

d(x 1 ,p 1 ,x 2 ,p 2 ,x 3 ,p 3 ) 



{X 1 ,P 1 }{X 2 ,P 2 }{X 3 ,P 3 } 



d(xi,p 1 ,x 2 ,p 2 ,x 3 ,p 3 ) 
{X U P 3 } {P U P 2 } {X 2 ,X 3 } - {Xt,P 2 } {A 2 , P} {X 3 , P 3 } - 
{X U P 3 } {A 2 , P 2 } {A 3 , P} - {Xi,Pi} {A 2 , P 3 } {A 3 , P 2 } + 
{Ai, A 2 } {Pi,P 3 } {A 3 , P 2 } + {Xx,P 3 } {A 2 , Pi} {A 3 , P 2 } - 
{A 1; A 2 } {P 2 , P 3 } {A 3 , Pi} + {A l5 P 2 } {A 2 , A 3 } {Pi, P 3 } - 
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{X U X 3 } {P 1; P 3 } {X 2 ,P 2 } + {X U X 3 } {X 2 , Fx} {P 2 , P 3 } + 
{Xt, X 3 } {Pi,P 2 } {X 2 ,P 3 } - {X U X 2 } {Pi, P 2 } {X 3 , P 3 } - 
{X u P x } {X 2 , X 3 } {P 2 , P 3 } + {Xx, P 2 } {X 2 , P 3 } {X 3 , P x } . (37) 

This expression can be checked with the help of the computer. 

It is easy to see that to obtain such a formula one needs to start 
from the term {Xx,P±} ■ ... ■ {Xjj,Pjj} and add to it all possible 
permutations. Factor multiplying each term is either +1 for even 
permutation and —1 for odd permutation. In the above expression all 
variables are ordered, i.e. in each Poisson bracket Xi is before Pj, Xi 
is before Xj if j > i, and Pi is before Pj if j > i. 

For deformation ([3HI]) Jacobian (f37|) simplifies to 

d(^Pi^ 2 ,P 2 ,X 3 ,P 3 ) &p22 2 

d(x 1 ,pi,x 2 ,p 2 ,x 3 ,p 3 ) 

In the case of small deviation of the deformed Poisson brackets 
(fTTj) from the canonical ones (/^ — 5^, hij, gij ~ 0) expression ([32]) 
can be simplified significantly in the linear approximation over these 
deviations. It is easy to see that in this case only the first term 
{Xi, Pi} • . . . • {Xd, Pd] contributes to the Jacobian and 

D D 

j = n mx, p) = i + p) - 1). 

i=l i=l 

Such an expression is useful for the analysis of low temperature be- 
havior. 
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Figure V. Temperature dependence of heat capacity of the harmonic oscil- 
lators ensemble per one particle. /? = /?' = 0.01, % = u = 2m = 1, T is 
given in hu units. Dashed line — heat capacity in the non-deformed case 
(/3 = (3' = 0), solid line — heat capacity calculated according to the direct 
definition of the partition function ([6]), dotted line — according to classical 
approximation (1T51) . 
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